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We develop a method to determine spatial configurations to realize decoherence-free subspaces 
for spatially extended multi-particle systems. We have assumed normal reservoir behavior including 
translational invariance of the reservoir and preparation in stationary states or mixture thereof and 
weak Markovian system-reservoir coupling that requires energy transfer. One important outcome of 
our method is a proof that there does not exist a multi- particle decoherence-free subspace in such 
systems except in the limit that the spatial extent of the system becomes infinitesimal. 

PACS numbers: 



I. INTRODUCTION 

The notion of a decoherence-free (or noiseless) sub- 
space (DFS) has been proposed as one strategy to com- 
bat the effects of decoherence in quantum computa- 
tion [1,0,0,01 and has the advantage of being passive, in 
contrast to quantum error correction, which requires ac- 
tive syndrome measurement and feed forward [l], 0, 0] . 
The construction of a DFS relies on the existence of 
symmetries in the decoherence process that allow certain 
subspaces of quantum states to be completely decoupled 
from the environment . A DFS is of practical interest be- 
cause it would reduce errors due to decoherence, thereby 
lowering the quantum error correction overhead in quan- 
tum information processing, and be especially useful in 
providing longer quantum memory. 

Extensive formal analyses of the DFS have been per- 
formed [I, 03, OH, [H, M, EI El and DFS have been 
applied in quantum information processing to reduce er- 
rors due to decoherence [l7], EH EH & UK- However, 
realizing a DFS (or approximate DFSs that would reduce 
decoherence over a finite, but long time scale) under gen- 
eral conditions is required. Our goal here is to address 
this need, namely to introduce general realistic principles 
for decoherence of systems of particles, to create a for- 
malism for describing these systems in order to determine 
systematically when a perfect DFS is impossible, and to 
demonstrate by example how an approximate DFS can 
be created even though a perfect DFS cannot be. 

We assume that the DFS is encoded into a collection 
of N identical particles, each with a D-dimcnsional in- 
ternal Hilbert space Ti. n spanned by the states 

{|a)„;a = l,. ..,£>}, n=l,...,N (1) 

Although the particles are identical, they arc not co- 
located (except when we refer to the "Dicke Limit" [22j 
in which the separation of particles is allowed to tend to 
zero). The position of particle n is denoted x„. In the 
absence of any coupling to the environment, each internal 
level \a) n is a stationary state of the system's Hamilto- 
nian evolution with eigenvalue u> a (we set K = 1). 



This system interacts with a reservoir, which is as- 
sumed to have no memory (Markovian). At time t = 0, 
the state of the system p§ and the state of the reser- 
voir are independent and then commence interacting. 
This separability of system and reservoir states is stan- 
dard in the theory of open quantum systems, and we 
assume that the evolution of the system in the presence 
of the reservoir can be described by a completely posi- 
tive map on the system, or equivalcntly by a Lindblad 
master equation (2g, H3 • We add the additional assump- 
tion that the reservoir is translationally invariant, which 
is reasonable but is also an important and restrictive as- 
sumption in studying spatially extended systems of parti- 
cles. The generator of translation is a vector momentum 
operator P. We assume that the Hamiltonian reservoir 
dynamics depend on |P| and not on the direction of P, 
i.e. the reservoir has no preferred direction. We refer to 
a reservoir satisfying these conditions as a homogeneous 
isotropic Markovian reservoir. 

Finally we place some reasonable restrictions on the 
interaction between system and reservoir. Specifically, 
we consider an interaction Hamiltonian of the form 

N 

H I = -Y,Y,<E i (xn), (2) 

n— 1 i 

with reservoir operators £^(x n ), which depend on the po- 
sition of the interacting particle n, and system operators 

W\ (3) 

that raise and lower the system energy levels. Here 
Lo a p = uj a — top ^ in the sum, which ensures that 
terms involving states |a) and |/3) of the same energy are 
excluded. Thus the reservoir absorbs and delivers energy 
to the system but does not contribute to strict dephas- 
ing. The Hamiltonian will be used to derive a master 
equation for the system. 

In summary, we have a system of particles with internal 
levels interacting linearly with a reservoir that is trans- 
lationally invariant, whose Hamiltonian is independent 



2 



of direction. These assumptions describe typical realistic 
systems, for example atomic or molecular gases, but have 
serious implications for constructing a DFS. On the one 
hand, we will find that we rule out the existence of a DFS 
under quite general conditions, but on the other hand, we 
will establish a foundation for obtaining a decoherence- 
supprcssed subspacc that is an approximate DFS. 



II. THE DECOHERENCE-FREE SUBSPACE 

The density operator for the system, traced over the 
reservoir, is given by p, which is a bounded operator on 
the Hilbert space Hs for the system of N particles, each 
with D internal levels. Under the assumptions given in 
Sec. HI we can express the evolution of p by the master 
equation 



with 



P = ~i[H e s,p] + L D [p], 



H eS = H s + A 



(4) 



(5) 



the Hermitian effective system Hamiltonian, where Hg 
corresponds to the Hermitian system Hamiltonian, A is 
the Hermitian contribution from the environment, and 
the nonunitary dynamics are incorporated into the deco- 
herence propagator 



the environment 0, H[ (i.e. p(t) = e ltHs pe ltHs or, equiv- 
alently p = — i[Hs,p] and —i[A,p] + Ljj[p] = 0). In 
Appendix [5] we show that for the systems we consider, 
this definition is a special restricted case of Definition [1] 

To determine necessary conditions for DFS we intro- 
duce 

Definition 2. An instantaneous pure decoherence-free 
subspace (IPDFS) Hipdfs is a subspace ofHs such that 
all density matrices p 6 D(Hipdfs) fulfill p 2 = p and 
L D [p}=0. 

Remark 4. We note that in general condition Lp)[p(t)] = 
at a specific time t does not guarantee unitary evolution 
of p(t) at a later time because H c g can potentially drive 
p(t) out of IPDFS. However, every state in DFS is an 
instantaneous pure decoherence-free state. Thus DFS is 
a subset of IPDFS. 

Now we establish a simple criterion for identifying 
whether a state is in an IPDFS and thus determine nec- 
essary condition for DFS, beginning with a definition of 
the decoherence operator. 

Definition 3. The decoherence operator for a system 
with decoherence propagator is 



M 

E 



XiJjJi. 



(8) 



1 



m 



Ld[p} = ^J2 x ^ 

1=1 



[Ji,pjh 



[Jip, # 



(0) 



In Eq. © , the real numbers A/ > are time-independent 
and operators {Ji}fi 1 form a subset of a complete basis 
for the space of bounded operators defined on Hs ■ 

We use the same definition for DFS as in j28|. Let 
D(Ti) be the set of all density matrices defined for a 
quantum system associated with the Hilbert space H. 

Definition 1. Let the time evolution of an open quan- 
tum system with Hilbert space Hs be governed by Eq. Q). 
Then a decoherence free subspace Hdfs is a subspace of 
Hs such that all density matrices p(t) 6 D(Hdfs) fulfill 



P = -i[H e s,p] \ft 



(7) 



Remark 1. Note that Definition [T] implies that Lo[p] = 
for all times t. 

Remark 2. The evolution of a state in 7Ydfs includes the 
environment-induced unitary evolution given by A, but 
not the decoherence effects. 

Remark 3. An alternative definition for a DFS requires 
that states experience unitary evolution due to the sys- 
tem Hamiltonian Hs only. In this case states are un- 
changed by (unitary and non-unitary) interactions with 



Remark 5. Note that the operator f is positive semidef- 

2 



> 



inite because V € Tts, (^P\^\% 
0. 

Inserting expression ([5]) for H c g into Eq. ^ and ex- 
pression ([5]) for the decoherence operator into Eq. ©, 
we see that general form of the master equation can be 
written as 

f - f 

p = - i[H s ,p] + (-iA - -)p + p{iA - -) 



(9) 



i=i 



if as the non- 



Henceforth we refer to the term iA 
Hermitian transition operator. 

The following proposition establishes necessary and 
sufficient conditions for the existence of an instantaneous 
pure decoherence-free state. A full proof is given in [28| . 

Proposition 4. For a pure state p = \ip)(t/j\, Ld[p] = 
iff Ji\ip} = ci\tp) and T\ip) = g\4'), where g = 

We demonstrate later that only operators J; that have 
no non-zero eigenvalues appear in the systems we con- 
sider in this paper. In this case the following proposition 
becomes very useful. 
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Proposition 5. Suppose that all eigenvalues of operators 
J i for I = 1, . . . , M are equal to zero; then the state 
is in an IPDFS iff it is an eigenvector of the decoherence 
operator T |3J) with an eigenvalue of zero. 

Proof. From Proposition U wc know that if \ip) is in 
an IPDFS then Ji\tp) = q \ip) for I = 1,...,M. The 
supposition for Proposition [5] implies that q = for 
1 = 1 M. Thus 



M 



1=1 



Conversely, if 



f|^) = o|v>, 



(with |||^)|| + 0), then 



M 



M 



= W\J2hJ}Ji\i>) = J2 x iP^ 



i=i 



(10) 



(ii) 



(12) 



Because A; > we can conclude that Ji\ip) — 0\tp) V I. 
Proposition 0] then guarantees that is in an IPDFS 

□ 

Remark 6. As the eigenvalues of the decoherence op- 
erator r are the inverse lifetimes of the corresponding 
eigenstates, Proposition [5] shows that a state l^) is in an 
IPDFS iff it has an infinite lifetime. 

One purpose of this work is to study whether or not 
an IPDFS (and consequently DFS) is possible for parti- 
cles located at different positions. We therefore have to 
distinguish IPDFS (DFS), which require two or more par- 
ticles from those that already exist for a single particle, 
i.e. from internal IPDFS (DFS) states. 

Definition 6. The single-particle IPDFS 'H^Ipdfs ( ~ 
Hn of particle n is spanned by all states which are 
decoherence- free if the system consists only of particle 
n, i.e., if the sum in the Hamiltonian (0) contains only 
one particle. The single-particle IPDFS for the complete 
N -particle system is the tensor product of the individual 
single-particle IPDFS, TL s -ipdfs = ®tl=i^ipdfs- 
Definition 7. The multi-particle DFS is the complement 

'Hm-IPDFS = 'H IPDFS \ 'Hs-IPDFS- 

The single-particle and multi-particle DFS can be de- 
fined similarly. The single-particle IPDFS (DFS) corre- 
sponds to a tensor product of IPDFSs (DFSs) for each 
particle. For instance, if all particles are atoms prepared 
in their ground states the corresponding iV-particle state 
is stable against decoherence due to spontaneous decay. 
The multi-particle IPDFS (DFS) typically involves states 
in which some or all particles are entangled and relies on 
destructive interference of the transition amplitudes for 
energy transfer from one particle to another. In the next 
section we introduce reasonable assumptions about real- 
istic reservoirs to study whether or not a multi-particle 
IPDFS (DFS) can exist. 



III. IPDFS IN REALISTIC SYSTEMS 
A. Markov-Born master equation 

In a realistic system, the locations and spatial config- 
uration of the particles are important. No two particles 
can occupy the same point, i.e. x„ ^= x„/ for n =/= n' . 
However, particles may be very close together and the 
limit that x„ = x n / V n, n' is referred to as the Dickc 
limit [22;]. The reservoir operators E l (pc) introduced 
above are position-dependent, and the reservoir Hamil- 
tonian TJr governs the dynamics of these operators. For 
example, the set {E 1 } corresponds to the vector compo- 
nents of the electric and magnetic fields for the case of 
an electromagnetic reservoir. 

The system+ reservoir interaction is given by Eq. 
which yields the master equation in the Markov-Born 
approximation: 

p{t)=-i[H s ,p]- T At ]T [[alM-^p 

n,m,i,j 

x (^(x„,r)^'(x m ,0)) + [pdU-^d. 
x (^'(x m ,0) J B i (x„,r) 



(13) 

In this master equation 

d\{-T) =cxp(-iH S T)d l n exp(iH S T) 

= £ a n ^e-^dl , (14) 

and E(t) = U^EUr are the operators in the interaction 
picture, with [/r = exp(— iH^r). 

B. Properties of the reservoir 

Master equation (fT"3"|) describes the evolution of N par- 
ticles interacting with a reservoir, whose propagation can 
be described by the correlation function 



GvL(t) =6(T)(E\x n ,T)E>(x m ,0) / 

=6(t)Tt r (^^(x„,r)^'(x m ,0)) , (15) 



with PR the density operator for the initial state of the 
reservoir. In what follows we assume that the reservoir 
has the following reasonable properties. The last two 
points below are standard assumptions in open system 
theory using Markovian master equations, which we state 
here explicitly for completeness. 

(i) is time-independent. 

(ii) The reservoir is initially in a stationary state: 
[pr, Hr] = 0. Then is diagonal in the -Hr- 
cigenbasis {|^ r k)}, with r and k quantum numbers 
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for these states (the meaning of k is explained in as- 
sumption (hi)) and {r2 r k} the corresponding eigen- 
values: 



Pr= d 3 k ^p r k|VVk)(VVk|, 



(16) 



with {p r k} the probabilistic weighting of these 
cigenstates. 

(iii) The reservoir is translationally invariant: for P the 
generator of translations for the reservoir in three 
dimensions. [i?R,P] = [pr,P] = 0, and the spec- 
trum of P is continuous. Typically P is the total 
momentum operator of the reservoir. The vector 
quantum number k thus has meaning as the mo- 
mentum: P|VVk) = k|-!/Vk)- 



(iv) .Hr is a function of 



and independent of the ori- 



entation of P. This assumption is widely valid. For 
example a reservoir of M interacting particles with 
Hamiltonian 



^-^ 2m 



(17) 



j>i 



has a corresponding generator of translations P = 
p;. Expressed in relative and center-of-mass co- 
ordinates, the Hamiltonian becomes 



2Mm 



H Te \, 



(18) 



with H TC \ depending only on the relative coordinates 
and not on the orientation of P. 

(v) Translations of the reservoir operators E i (x) are 
generated by the operator P, 

£ l (x + x) = c lp - x E l (x a )e-* p -* . (19) 

(vi) For technical reasons we need to introduce the set 
of k vectors 



K Ul r lS (k') = {|k| : fl rk = n s 



■w)} 



(20) 



and require that the reservoir transformation func- 
tion 







.sk' 



E l (0) 



Ei(0) 



sk' 



|k|eiW,s(k') 



(21) 



is either (a) a polynomial in terms of Cartesian com- 
ponents of k, or (b) is equal to f(k)A l J s such that 
A l J s independent of k and the support of f(k) con- 
tains |k| = r for some r £ K Ut r iS (k'). The transfor- 
mation function describes the change in the state of 
the reservoir when a transition in the system occurs 



that changes its energy by h\o. The expression in 
Eq. (f2"Tj) is constrained to the set K^^^lk 1 ) because 
only transitions that conserve energy in the overall 
system (quantum system+reservoir) are important. 
Case (b) describes a reservoir composed from scalar 
particles. Note that there exists k G K^^^ik') for 
some oj, r, s, k' such that the transformation func- 
tion (|2"Tj) is not identically equal to zero for |k| = k, 
since otherwise the quantum system would not be 
coupled to the reservoir. 

(vii) All energy-nonconserving transitions in the system 
of N identical particles average to zero on the time 
scale relevant to the Markovian master equation. 
This condition guarantees that the Markovian mas- 
ter equation we derive is completely positive (i.e. 
physical) j29| . Inclusion of energy-nonconserving 
transitions modifies A [13, HH, HH , which is irrele- 
vant for the purpose of this paper since, as we have 
shown, only the decoherence operator T is needed 
for determining IPDFS. 

(viii) System+reservoir coupling is sufficiently weak so 
that second-order perturbation theory is valid. 

(ix) The reservoir correlation function decays rapidly 
in time so that backreaction is effectively instanta- 
neous, which underpins the Markovian master equa- 
tion and guarantees that the frequency spectrum of 
G contains only frequencies that are much larger 
than the coupling parameters. Consistency of the 
master equation then implies that the processes un- 
derlying it are energy conserving [29| . 

Definition 8. A homogeneous isotropic Markovian 
reservoir is any reservoir that satisfies all conditions (i) 
to (ix) above. 

From these assumptions, the Fourier transform of the 
correlation function G is 



g&»=e / d3k ' / d3k ^k<c 



i(k-k')-(x„-x ra ) 



Vvk) (Vv 



E*(0) 



x I ttS(uj — f2 r ksk') + i- 



Ei(Q) 
V 



4's 



(22) 



CO — fi r ksk' 

Here 

^rksk' = ^rk — ^sk' (23) 

and details of the calculation are provided in AppcndixlBl 

C. Decoherence operator and correlation function 

After comparing Eqs. @ to (fTS")) and recalling the defi- 
nition for the correlation function GJf m (r) from Eq. (|15|) . 
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one can identify the non-Hcrmitian transition operator 
specific to our system as 



Eq. (|25|) and retaining only energy-conserving terms, we 
obtain 



^ + ¥ = r drX:^ ro (r)44(-r). (24) * A + ^ = E/d 3 k' /d 3 k £ T ra (k,k> a/3 ) 



From Eq. (|14|) . we can rewrite Eq. (j24[) as 



drGl(r)e- ! ™^'. 



(25) 



By comparing the integral in Eq. (|25p to the integral in 
Eq. l|B2p . one sees that 

dr Gjf m (r)c-^'»' = Gjf m (-o;^ ) . (26) 

Now we introduce the collective operator 
S rs (k, k', w) = X e i ( k - k ')' x »E rs; n (k, k', w) , (27) 

with 

E rS ;n (k, k',w) = X &n,aPy/Psk' 

a, P 

X (VVk £i<(0)d* a/3 



V^k' ) ■ (28) 



Here S rs;n (k, k',w) is an operator that excites (w > 0) 
or de-excites (u> < 0) the n th particle. 

Operator (|2T|) is a direct sum of operators on each par- 
ticle and incorporates the reservoir state and coupling 
term to create a specific raising or lowering operator. 
The operator is defined on states \a) and \/3) connected 
by raising (or, if lu < 0, lowering) operators a nt0t p so that 
the energy difference between these states is lu. Some- 
what analogous to the number operator, we combine S 
with its adjoint to define 

f „, „ (k, k', u) = H s . „ (k, k', lu) ± rs , „ (k, k', lu) (29) 

for the n th particle, and 

T ra (k, k', u) = ±U (k, k', w) t rs (k, k', lu) (30) 

globally. 

Henceforth, we impose assumption (vii) of Sec. IIIIB1 
which in this context requires that all terms in Eq. (|25|) 
that fail condition Lu a p = to a i pi can be neglected. This as- 
sumption ensures the consistency of the Markovian mas- 
ter equation [2{|. Then by combining Eq. (|2"2"|) with 



x n6(ijj a p - Orksk') + i- 



aft ^rksk' 



(31) 



The decoherence operator Y is the Hermitian part 
of the right-hand side of Eq. (|3"Tj) . Because operator 
T rs (k, k', LO a p) is Hermitian by construction, it follows 
that the term involving the Dirac <5 generates the Her- 
mitian decoherence operator, whereas the term involving 
the principal value V generates the energy level shifts 
(the "Lamb shift" in quantum optics). We thus find 

t u =2tt^ J d 3 k' j d 3 kT rs (k,k» 

r.s 

x5(cj-n rk3 v). (32) 



and 



(33) 



By assumption (iv) of Sec. lIIIBl the frequency difference 
^rksk' depends on the absolute value of k only. Therefore 
the Dirac 6 fixes the absolute value of |k| by the condition 



n 



rksk' 



(34) 



that is, the energy differences in the reservoir Q r k — O s k' 
equal the corresponding energy differences for the system 
transition uj a — top. For a given vector k' and energy 
difference cu a p Eq. (|34|) determines a set of values |k| for 
which it is fulfilled. We denote this set by 



A^(k') = {|k| : lU = fi sk , +«;)}. 



(35) 



which is not empty for some tv; otherwise the system and 
reservoir would never exchange energy. 

Let iSfc denote a sphere of radius k. Then the integral 
with respect to k in the expression for the decoherence 
operator T is reduced to an integral over the surface of the 
sphere Sk with k 6 K^^^Ck'), where lu takes all possible 
frequency differences within one particle: 



lu). 



r = 2lT E / d3k ' E / d2k ( k ' k ' 

(36) 

Although the domain for quantum numbers r, s, and k' 
could be discrete, continuous, or a combination thereof, 
assumptions (i-ix) imply that k € Sk, which is a sphere of 
radius k. We have thereby shown that decoherence oper- 
ator L is related to T and determined the general expres- 
sion for S. This information will allow us to demonstrate 
in the next section that multi-particle IPDFS (and conse- 
quently multi-particle DFS) cannot exist in this system. 
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IV. CONDITIONS FOR MULTI-PARTICLE DFS 
STATES 

We seek a necessary criterion for the existence of multi- 
particle DFS states. We do this by establishing necessary 
and sufficient conditions for existence of multi-particle 
IPDFS for spatially separated particles in a homogeneous 
isotropic Markovian reservoir. We start with excluding 
the single-particle IPDFSs and we prove the following 
Proposition. 

Proposition 9. There is no multi-particle DFS for 
spatially separated particles in a homogeneous isotropic 
Markovian reservoir. 

Proof. In this proof wc make use of Proposition^ There- 
fore we need to check whether our system satisfies the 
hypothesis of Proposition [5j i.e. whether all eigenvalues 
of Ji are zero. Each operator r w can be decomposed as 
in Eq. ([5]). Index I = (n, to, a, (3) for operator J; de- 
pends on n (particle index), on oo (possible frequency 
differences within one particle), and on energy levels la- 
beled by a and /? such that co a p is constrained to equal 
oo. Each operator J\ is some linear combination of opera- 
tors a n ^ a p with a and (3 constrained so that oo a p is some 
fixed non-zero value. This means that there exists a rep- 
resentation with all operators Ji either upper or lower 
triangular matrices with zeros along the diagonal. Then 
we know that all eigenvalues of operators Ji are equal to 
zero and Proposition [5] applies to our system. 

First we need to determine the condition for existence 
of single-particle IPDFS states and remove these from 
the subspace Hipdfs ■ According to the definition [7] and 
Proposition [5l we are seeking states \a) € H n such 

= 0, with r„ the decoherence operator T 

particle. At the same time, we know 



that 



restricted to the n 
that 



f „ = 2tt Y, I d 3 k' Y i d2kt -; " ( k ' k '' u ) 



As a result, a state |a) G ?s™pdfs = ^« n ^ipdfs iff 



S rs; „ (k,k', u>) | a) 







(37) 



for all the external quantum numbers r, s, k', for all 
possible frequency differences w within one particle, for 
all values k in the set K^^sO*-'), and for all vectors k that 
lie in the sphere Sk- The space of single-particle IPDFS 

of the total system is then H s -ipdfs = ©^i^ipdfs 

Now we determine under what conditions the multi- 
particle IPDFS W m -IPDFS = ^ipdfs \ H s -ipdfs is not 
empty By an analogous argument we infer that 



la) e Hipdfs 



t r , (k, k',co)\a) 



= 



(38) 



for all quantum numbers r, s, oo, k', k and for all vectors 
k that are in the sphere Sk- 



Let \yi), ■ ■ ■ , \ yN>) be a complete basis for 7i m -iPDFS- 
Then 

Ye i ^- k '>^E rs ., n (k,k',u>)\y j )=0 

n 

Vi = l,...,7V', (39) 

with TV' = dim(7i m _iPDFs)- By assumption not all S rs: „ 
do annihilate the state \yj) because it then would be an 
element of 7Y s -ipdfs- 

Eq. (f3"9")l implies in particular that 

V^ e i(k-k')-(x„-x m ) 



X ( yj \±l m (k,k>,Lj)£ rs;n (k,k',U,)\ yj ) 



|k|ei^, r ,»(k') 



(40) 



for all j = 1, . . . , TV', for all quantum numbers r, s, oo, 
k', k, for all vectors k that are in the sphere Sk, and 
for all m = 1,...,N. Note that Eq. |gQJ) is the ker- 
nel T rs (k, k',u;) needed to calculate operator T as in 
Eq. (f36|) . Condition (j40|) is not trivially satisfied because 
for each j there is at least one m such that 



Y rs . m (k,k', oo) 



t rs . m (k, k',oo) \yj) 



since \yj) £ H s _ipdfs- Let 



^0, 
(41) 



^rs; m s - j k', Oj) S rs; n (k, k'. 



00 



II j 



(42) 



We note that there exists k G K^^^ik') such that c™ n ^ ^ 
since otherwise quantum system would not be coupled 
to the reservoir. Then 7i m _ipDFS is non-empty only if 







(43) 



is satisfied for all possible vectors k that are in the sphere 
Sk, for all j = 1, . . . , TV', and for all m = 1, . . . , TV with 
not all j =0. Now we make use of the two alternative 
assumptions (vi, a) and (vi, b) in Subsec. IIII Bl 

Assumption (vi, a):- This assumption allows us to con- 
clude that 



c" ■ = h n 



(k)&: 



m,ji 



(44) 



with j (k) a non-zero polynomial function in compo- 
nents of vector k defined on the sphere Sk and j co ~ 
efficients independent of k. Then Eq. (|4"3")) becomes 



(45) 



which is equivalent to determining conditions when a set 
of functions 



{h» 7 (k)c ikx "}£Li 



(46) 
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defined on a sphere Sk can be made linearly dependent. 
It is well known that functions of this form are linearly 
independent unless there exists m 7^ n such that x„ = x m 
(see Appendix [C")l . Thus requirement (|4T))) can only be 
fulfilled if 3 m 7^ n such that x n = x m , i.e. at least two 
particles occupy the same point in space. 
Assumption (vi, b):- This assumption yields 



/(k)C 



(47) 



and {bm,j} are coefficients independent of k. 

As the function /(k) is not zero for k G iSfc, condition 
(|43|) is equivalent to requiring that 







(48) 



for all possible vectors k <G Sk with not all j = 0. This 
is a special case of the situation we considered above and 
thus we reach the same conclusion as before: a multi- 
particle IPDFS can exist only when at least two particles 
occupy the same point in space. Consequently, a multi- 
particle DFS also does not exist when no two particles 
occupy the same point in space. □ 

Although Proposition[5]rulcs out a multi-particle DFS, 
states that do not decohere over a long lifetime, i.e. sig- 
nificantly enhanced beyond the longest lifetimes of the 
individual particles, may be just as practical. Because 
of the nature of the reservoir being continously transla- 
tionally invariant (assumption (iii)), we expect that this 
significant enhancement is achieved when particles are 
close together on some scale determined by the reservoir; 
we explore this case in the next section. 



V. DECOHERENCE SUPPRESSION 

We now apply the formalism developed in previous sec- 
tions to a well-studied model in quantum optics, namely 
a set of two-level atoms interacting with a multi-mode 
radiation field, which serves as a basis for studying col- 
lective phenomena such as superradiance and subradi- 
ance [22|, [23|, \Z% [25| . In fact some of the notation was 
inspired by this example, such as using the variable d, 
which is typically used for the electric dipole moment 
of the atoms, and E l , which is used for components of 
the electric field. In general the reservoir for a collec- 
tion of atoms, such as a gas, closely satisfies the ho- 
mogeneous isotropic Markovian conditions, and the sys- 
tem+reservoir coupling is as described in previous sec- 
tions. 

This model also provides a welcome simplicity in that 
each atom contains only one excited state so an internal 
DFS is not possible; hence there is no non-trivial single- 
particle DFS (i.e. other than the ground state for the 
system). As atoms cannot be truly co-located, we know 
from Proposition [5] that in this system a perfect DFS 
exists only in the Dickc limit for which atoms are co- 
located. 



The question we address now is whether an approx- 
imate DFS can exist, and how well it behaves, in con- 
ditions that are close to the ideal condition required for 
the DFS to exist. The atomic case gives us a concrete 
and well-studied example that gives insight to answer this 
question. 

Let us begin by assuming that the two-level atoms are 
located at positions Xj (j = 1, . . . , N) and that atomic 
dipole moments are given by d. Then the decoherence 
operator has the form [H, [3J, HH, [H, [37j 



N 



(49) 



with 




sin k^x 



k Xjk 

coskoXjk sin fco fc 
(k Xjk) 2 (koXjk 




+ 



and k = 27r/Ao, x.,- fc = Xj - x fe , x jk = \x. jk \, a k is the 
lowering operator for the k th atom, 70 is the Einstein A 
coefficient, and Ao is the resonant wavelength. We also 
define reduced decoherence matrix to be *7* = (7jfc). 

Eigenvalues of the decoherence operator, T, provide 
inverse lifetimes of the states in the system. For many 
applications, it is sufficient to look at the reduced matrix 
< 7 > . This matrix describes the states in the system of N 
two-level atoms with one atom in the excited state. It 
also describes the states for N two-level atoms with all 
atoms except one in the excited state (see Appendix [D]) . 
In general, it is much easier to deal with an N x N matrix 
rather than with an 2 N x 2 N matrix representation 
of the decoherence operator T. 

For special cases where the location of atoms can be 
described with one variable, one can produce two dimen- 
sional plots that show the dependence of lifetimes on the 
separation between atoms. As an example, consider N 
atoms in a line. The separation between adjacent atoms 
is taken to be equal (a^i+i = r). Figs. l(a, c, d) show 
particular configurations of two-level atoms in a line as 
well as the orientation of the dipole moment for the atoms 
that we consider here. Figs. 2(a, c, d) depict the inverse 
lifetimes (eigenvalues of matrix ^7*) of the states in the 
system of two and four atoms, when only one atom is 
excited, as a function of r. In Fig. 1(b) we consider the 
system of four atoms that form a square. The inverse life- 
times for states that arise when only one atom is excited 
in this system are shown in Fig. 2(b). From the plots one 
can see that there are long-lived states (i.e., inverse life- 
times close to zero) when adjacent atoms are separated 
by less than a quarter of the emission wavelength (Ao/4). 
For each cluster of closely located atoms, one could use 
the states with longest lifetimes (smallest eigenvalues) to 
encode a qubit. 
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So how much advantage can one gain if the separation 
between atoms is small? In Fig. 2(c) the longest lifetime 
possible is 109 times longer than the lifetime of the state 
of an isolated atom when the separation between adja- 
cent atoms is just Aq/4. However, the second longest 
lifetime possible in this system is just 4.5 times longer 
than the lifetime of the state of an isolated atom. Thus 
if we have chosen to construct a qubit from these two 
states, our effective advantage would be just the lifetime 
of the second state. The situation is a little better if we 
consider configurations that have more symmetry: for 
example, consider the square in Fig. 2(b). Here the life- 
times for the two collective states are comparable when 
the separation between atoms is small. When the side of 
the square is equal to Ao/4, the two largest lifetimes are 
4.6 and 5.1 times longer than the lifetime of the state of 
an isolated atom. 
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FIG. 2: Inverse lifetimes (relative to 70) for configurations 
shown in Fig. 1. The two plots in (a) correspond to the 
configuration of two atoms in fig. 1(a). The dash-dot line 
corresponds to the symmetric state and the solid line to the 
antisymmetric state. The plots in (b)-(d) correspond to the 
configurations in fig. 1(b)- (d) respectively and the four lines 
correspond to the collective states that arise when one atom 
is in the excited state and all other atoms are in the ground 
state. 



FIG. 1: (a) The system of two two-level atoms in a line 
with each dipole moment co-aligned with the axis joining the 
atoms, (b) The system of four two-level atoms that lie on the 
corners of the square with each dipole moment co-aligned with 
one specific side of the square, (c) The system of four two- 
level atoms in a line with each dipole moment co-aligned with 
the axis connecting the atoms (d) The system of four two- 
level atoms in a line with each dipole moment perpendicular 
to the axis connecting the atoms. 

We have shown how lifetimes for quantum states in- 
crease when the separation between atoms is small for a 
collection of a few atoms. Similar effect can be observed 
in a large collection of atoms. Spontaneous decay rate 
from a collection of a large number of two-level atoms 
was recently studied in two configurations: (a) N atoms 
in a linear chain with the same distance between neigh- 
boring atoms [4ll |42| and (b) a circular configuration, 
in which atoms were placed equidistantly on the circum- 
ference of the circle [4J| ■ 

For the case of the linear chain, a numerical study 
showed that for large number of atoms (40-100) [4l|, sig- 
nificant decay rate suppression was observed when neigh- 



boring atoms were separated by a distance d < Xq/2. In 
the limit of TV — * 00, the decay rate was seen to have 
a jump-like behavior at Xq/2. This observation implied 
that the spontaneous decay rate is unstable and suscepti- 
ble to environmental perturbations in the vicinity of this 
critical value. 

Analysis of the circular configuration of atoms offered 
more surprising results. It was observed that for a given 
radius r of the circular configuration, once a certain crit- 
ical number of atoms is reached the lifetime of the maxi- 
mally subradiant state increases exponentially and is now 
stable to perturbations. For example, for r = Xq, the 
critical number of atoms is fourteen, and, at this mo- 
ment, the smallest interatomic distance between atoms 
is 0.45 x A . 



VI. DISCUSSION 

Proposition [5] states under which conditions a DFS 
does not exist, namely whenever the set of atoms is away 
from the Dickc limit. To find a DFS outside of the 
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Dicke limit we therefore have to find systems that vio- 
late these conditions. One way would be to consider a 
single-particle DFS with many such particles in the sys- 
tem. Another way is to engineer a reservoir that is not 
of homogeneous isotropic Markovian type. 

Although the single-particle DFS is excluded for two- 
level atoms, a DFS can be created by a collection of par- 
ticles that have a nontrivial internal structure. Unfortu- 
nately the spherical symmetry of atoms induces selection 
rules for angular momentum eigenstates that prohibit the 
existence of a DFS in isolated atoms. However, instead of 
atoms one can use molecules that do not posses a spheri- 
cal symmetry. A DFS in molecules is known in the liter- 
ature as spontaneous emission cancellation. It has been 
reported for Rydberg states in Na2 [3|| but the observa- 
tion has been questioned [iol ]. 

The conclusions we have reached about the existence 
of DFS in this work depend in a fundamental way on the 
fact that the quantum system was embedded in three- 
dimensional space. Requirements for existence of DFS in 
one-dimensional structures are different. The logic be- 
hind Proposition [5] tells us that in this case the sepa- 
ration between atoms is no longer required to be zero 
for the existence of a DFS. (In one-dimensional struc- 
tures, the reservoir transformation functions are reduced 
to scalars and condition for linear dependence is always 
satisfied.) Therefore, a multi-particle DFS may exist for 
atoms placed in an effectively one-dimensional structure 
such as a waveguide or an optical fiber. P. Zanardi and 
F. Rossi have discussed similar results in the context of 
semiconductor structures [HI]. 

In summary, we have developed a new way to de- 
termine the existence of both single-particle and multi- 
particle DFS, for realistic systems described by a Marko- 
vian master equation. We have demonstrated that a 
perfect multi-particle DFS requires co-located particles 
(i.e., the Dicke limit) for systems placed in a homoge- 
neous isotropic Markovian reservoir, which is not possi- 
ble. Our analysis shows however that it is nevertheless 
possible to have a single-particle DFS involving many 
atoms. Also a multi-particle DFS for atoms situated in 
one-dimensional structures might be possible. Further- 
more we have established here a foundation for studying 
approximate DFS and applied it to a set of N two-level 
atoms interacting with a multi-mode radiation field. 
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APPENDIX A: COMPARING TWO 
DEFINITIONS FOR DFS 

Here we show under which condition different condi- 
tions for DFS arc equivalent for systems coupled to a 
homogeneous isotropic Markovian reservoir. 

Proposition 10. For a pure state p = {tp\ of a system 
coupled to a homogeneous isotropic Markovian reservoir, 
if-i[A,p]+L D [p] = then L D [p] = 0. 

Proof. Suppose that — i[A, p] + Lo[p] = holds. We then 
can evaluate the expression 



= M - i[A,p]+L D \p]\1>) = ^\L D [p]\i>). (Al) 



But 



M 



\Ld[ p ] iv>> = EM (v^w 



7 = 1 



(A2) 

Ji\ip) can generally be written as Ji\i/j) = ci\tp) + 
with \ipj~~) some (non- normalized) state that is orthogonal 
to the state If we substitute this into Eq. (|A2[) we get 
E^i M(i>i-\i>i-) = 0. Because A; > for I = 1, . . . , M 
wc find \\\ipA)i\\ = 0, i.e., J,|^) = ci\i)) : Thus |t/>) is 
an eigenstate of all the error generators J; appearing in 
Eq. [5] As we have argued before for the systems im- 
mersed in the homogeneous isotropic Markovian reser- 
voir, operators J/ can have only zero eigenvalues. Using 
Proposition 0] we can infer that Lr>[p] =0. □ 



APPENDIX B: FOURIER TRANSFORM OF THE 
CORRELATION FUNCTION 

We derive Eq. (|22|) starting with Eq. (fTS"]) . Expanding 
G^ m (r) in terms of |^Vk) an d introducing the unitary 
operators C/ x „ = exp(ix„ • P) results in 



GL(r) =0(t)J2 J d 3 k'| d 3 k p sk , 



X (lp r k |#?(x m ,0) -0 sk . 

=9{t) I d 3 k f d 3 k' p s v 

r,s 

x (V> s k< \ululE l (0)U^U R 
x (l>*\ul n E i (P)U Xn \l>*' 

^( T )E / d3k '/ d3k Psk '° 







rk 



z(k— k') ■ (x„ — x„ 



-ii-^r-ksk' 



x U rk Ei(0) V> 8 k 



(Bl) 
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We want to calculate the expression for 

Gt» = f°° drG^ n (r)cxp( lW r). (B2) 

J — oo 

The time-dependent part of the Green's function G^ m (r) 
in Eq. (|Blj) is of the form e _?Tn ' ksk ' . Thus the integral 
with respect to r in the expression for G^ m (ui) is 

/>00 

lim / drc' T(w - !! ^' ) -»' T 
v-° Jo 



»7— >o uj — fi r ksk' + iv 
= (tt5(uj - n rksW ) + i £ ) . (B3) 

\ UJ — H r kaW J 

Then the expression for the Fourier transform of the 
Green's function is 

X (^, k , |i*(0)|Vrk> (^k|^(0)|^k'> 

x ( vr(5(a; - fi rks k') +t pr )■ (B4) 

\ U — Sirksk'/ 

We use this expression in Eq. (|22p. 

APPENDIX C: LINEAR INDEPENDENCE 

Here we verify the claim made in the proof of Proposi- 
tion [9] that the set of functions {Pj(k)e Xn }j—i, where 
Pj (k) are polynomials in Cartesian components of k de- 
fined on the sphere of radius r (|k| = r), are linearly 
independent if x„ ^ x m for n ^ m. 

Before we proceed, we need to introduce some nota- 
tion. Let components of the vector k be denoted by 
k\, fc 2 , k% and let a = (ai, a 2 , CK3) be a collection of non- 
negative integers. We define k Q = fc" 1 /^ 2 ^ 3 . Then a 
polynomial P(k) is denoted as P(k) = J2 a a ak a , where 
coefficients a a are non-zero only for finite collection of a. 
Let \a\ = ct\ + 0L2 + Qi3- The /it?/ degree of the polynomial 
P(k), denoted by dcg(P), is the largest \a\ of all terms 
in P(k) with non-zero coefficients a a . 

Proposition 11. Let Pi(k), P 2 (k), • • ■ , Pv(k) 
be non-zero polynomials in Cartesian compo- 
nents of vector k. T/ien t/ie set 0/ functions 
{Pi(k)e ?;kxi , P 2 (k)e ikX2 , • • • , Pv(k)e ikXN } is lin- 
early independent when vectors with real entries 
Xi , X2 , • • • , xjv are distinct. 

Proof. Functions Pj(k)e lkx ™ j = 1,...,N are linearly 
independent, if they cannot be expressed in the form 

ciPi(k)e ik - Xl + c 2 P 2 (k)e lkX2 +••• 

+c w Pv(k)e lk - x " = (CI) 



with Cj (j = 1, . . . , N) constants not all equal to zero. 

We assume that Eq. (|C1|) holds. Then we use mathe- 
matical induction to show that all Cj's are 0. 

For N = 1, we have ciPi(k)e ik ' Xl = 0. Since 
Pi(k)e lk ' Xl is a non-zero function, c\ = 0. 

Now suppose the statement is true for N functions, i.e. 

if 

ciPi(k)e ikxi + c 2 P 2 (k)e lkX2 + ■■■ 

+ Cw P A r(k)e lk - x « = 0, (C2) 

then Ci = for 1 < % < N. 
We want to show that if 

ciPi(k)e ikxi + c 2 P 2 (k)e lkX2 + ■■■ 

+c w+1 P A r +1 (k) e lk - x «+ 1 =0, (C3) 

then Cj = for 1 < % < N + 1. 

Let x„ = (x„ ; i, x„ ;2 , x„, 3 ) and let m n = deg(P„) + 
1. Consider the linear differential operator L n = 
Sj=i(gf" — * x n,jl) m " an d observe that 

L n (P„(k)e*' x ») = 0. (C4) 

For I + n, L n (P ; (k)e« kx ') = Q Ml (k)e ik x ' , (C5) 

where Qz,n(k) is a non-zero polynomial with deg(P) = 
deg(Q ;> „). 

The statement (|C4|) holds because the power of the 
term ( — ix nj l) in L n exceeds the degree of variables 
k\, fc 2 , &3 in polynomial P n (k) at least by 1. The state- 
ment (|C5p is verified later. 

We apply a linear operator Ljv+i to the Eq. (|C3|) . 
Then 

ciQi,AT + i(k)e jk - Xl + c 2 Q 2iJ v + i(k) e ik - X2 + • • • 

+c N Q N . N+1 (k)e lk ^ = 0. (C6) 

Polynomial Q;,jv + i(k) is non-zero and it has the same 
degree as P;(k) for 1 < I < N. Therefore by in- 
duction hypothesis c\ = c 2 = ■ • • = cn = 0. Thus 
CAr + iP/v + i(k)e ik ' Xjv+1 = 0, which is only possible when 
cn+1 = 0. 

In order to check that statement (|C5|) is correct, we 
introduce polynomials Qj (k) such that 

jjjT '^...'j fl(k)e*' x '. (C7) 

Notice that Q/,„(k) = J2j=i Qj, n ( k )- Sincc X; ^ x ™' 
there exists j G {1, 2, 3} so that x;j ^ x nj . We claim 
that deg(Q; n ) = deg(Pj). Indeed for m n = 1, 
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Therefore Qj (k) has the same total degree as P; (k) . 
Suppose the claim is correct for m n = M. Then for 
m„ = M + l, 

x ^(k)+i(xy-x n ,,)P i (k)) e ik x '. 

(C9) 

By assumption deg(Q[ n (k)) = deg(ffi(k) + i(x tij - 
x„ J )P i (k)) and deg(§g(k) + -x^P^k)) = 

deg(P(k)). 

Also if polynomial p(k) is constant, i.e. P(k) = 
c 7^ 0, then polynomial Qz in (k) has the form Qz,n(k) = 
c Ej=i (*( x *,j - x «.j)) m "- Since x; 7^ x„ for Z 7^ n, poly- 
nomial Qi ,„(k) is also non-zero. □ 

APPENDIX D: RELATIONSHIP BETWEEN 
DECOHERENCE OPERATOR AND REDUCED 
MATRIX. 

We analyze here the relationship between the reduced 
matrix and the decoherence operator T. In the Propo- 
sitions below we demonstrate the relationship between 
the eigenvalues of the matrix 7 and the decoherence 
matrix T, as well as establish the connection between 
eigenvectors for these two matrices. 

Proposition 12. // A is an eigenvalue of < t* ; then 70A 
is an eigenvalue of decoherence matrix T. 

Proof. Let (x\ ■ ■ -xn) 1 be an eigenvector of matrix *7*. 
Then 7 (xi ■ ■ ■ xn) = A (x\ ■ ■ ■ xn) for some real num- 
ber A. In terms of components this can be written as 

Sfe=i7ife a; fe = Xxj. 

We would like to show that 70 A is an eigenvalue for 
operator T. This will be true iff there exists state such 
that (V>|(r/7o — A/)|t/>) = 0. Consider a normalized state 
= x j|1>0j where the state \l,l) corresponds to 

I th atom in the excited state, while all other atoms are in 
their ground state. The state |0) denotes the collective 
ground state for all N atoms. Then ak\ip) = SmXi\0) and 
(ip\ajd-k\ip) = XjXk- Thus 

N 

(V>|(f/ 7o - MM =(VI( E wj<W) - A 

N 

= E 7jfe(V'|o']o'fc|V') - A 

= E 'YjkXjX k - A = (Dl) 



□ 

Proposition 13. Eigenvalues of the matrix 70^ are 
inverse lifetimes for collective states of the system with 
N — 1 two-level atoms in the ground state and one atom 
in the excited state. 



Proof. Lifetimes for N asymptotically separated atoms 
is the same as for N isolated atoms. Thus inverse life- 
times for collective states with one excited atom (all other 
atoms in the ground state) are equal to 70. There will be 
N such states. Matrix ^7* has N eigenvalues and for infi- 
nite separation between atoms its eigenvalues are 1. □ 



Proposition 14. If X is an eigenvalue of^y*, then (N — 
2 + A)7o is an eigenvalue of the decoherence operator T 
in the system of N two-level atoms. 

Proof. We proceed as in the proof for Proposition[T2"l Let 

[x\ ■ ■ ■ xn) T be such that X^fcLi ljk%k = Xxj. 

Consider the normalized state \ip) = Y^iLi a^z |0, Z) , 
where state |0, I) corresponds to I th atom in the ground 
state, while all other atoms are in their excited state. 
The state |1) denotes the collective excited state for all 
N atoms. Then &l\ip) = 5uxi\l), (ip\&jak\ip) = x.jX k for 
k ^ j and (V>|<rj>s#) = £^ fc=1 x ) = 1 ~ x l- Thus 

M (f /7b - (N - 2 + A)/) |V) 
= £ ljk (iP\a]* k \i>)-(N-2 + \) 

j,k=l 

N N 

= E ljkXjX k +N-2j24-(N-2 + X)=0 (D2) 

j,k=l k=l 

□ 



We have shown how knowledge of the reduced decoher- 
ence matrix ^ allows us to determine 2N eigenvalues of 
the full decoherence operator T. We also were able to 
derive the expressions for the corresponding eigenstates, 

which are \ip n ) = £j =1 tf\0,l) an d \(j> n ) = Eili ^"l 1 ' 
with n = 1, . . . , N, where |1, 1) is the state with only / th 
atom in excited state, |0, 1) is the state with only I th atom 
in the ground state, and (x™ . . . x r ^) T is the n th eigenstate 
of matrix 7 . 
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